A Re-examination of Quenches in 4 He 
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In the light of recent difficulties in observing vortices in quenches of liquid 4 He to its superfluid 
state we re-examine the Zurek scenario for their production. We find that the standard prediction 
is suspect because of a confusion over correlation lengths. 
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In thermal equilibrium the behaviour of simple systems 
experiencing a continuous phase transition is generic, as 
manifest in the utility of Landau-Ginzburg theory. Is this 
also true dynamically? The relevance of this question is 
that the early universe proceeded through a sequence of 
phase transitions whose consequences are directly observ- 
able, but whose detailed dynamics is unknown. Although 
it is difficult to measure an order parameter as it changes, 
many transitions generate topological charge or topolog- 
ical defects which can be detected. Motivated in part by 
Kibble's mechanism for the formation of cosmic strings 
in the early universe [0J^], Zurek suggested j|] that we 
measure the density of vortices produced during a pres- 
sure quench of liquid 4 He into its superfluid state, as well 
as the variance of superflow velocity. 

The scenario, as proposed by Zurek, is very simple. It 
is exemplified by assuming that the dynamics of the tran- 
sition can be derived from an explicitly time-dependent 
Landau-Ginzberg free energy of the form 



F(t) = / d 6 x 
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In (|I|) (j) — (0i + *02)/ V% is the complex order-parameter 
field, whose magnitude determines the superfluid den- 
sity. We identify a(t) as an externally driven time- 
dependent chemical potential. In equilibrium at temper- 
ature T, in a mean field approximation a(T) takes the 
form a(T) — ctae(T c ), where e = (T/T c — 1) measures 
the critical temperature T c relative to T. In a pressure 
quench at approximately constant T, T c will vary with 
time t, and we assume that e can be written as 



e(t) =e - —0(t) 
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for — oo < t < tq(1 + eo), after which e(t) = —1. eo = 
(T/T™ — 1) measures the original critical temperature 
T* n against the temperature T at which the quench takes 
place, and tq defines the quench rate. The quench begins 
at time t = and the transition from the normal to the 
superfluid phase begins at time t — €qtq. 

With £q = h 2 /2mao and tq = h/ao setting the funda- 
mental distance and time scales, the equilibrium correla- 
tion length £,(At) and the relaxation time r(At) diverge 
at the relative time At = t — €qTq = as 



/AtV 1/2 

f (A*) = eo^J , 
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As we approach the transition, eventually the relaxation 
time will be so long that the system will not be able 
to keep up with the temperature change. We estimate 
the time tz (and the relative time Atz = tz — £o t q) 
at which the change from equilibrium to non-equilibrium 
behaviour occurs by identifying r(Atz) with — Atz i-S. 
—Atz — \/toTQ . After this time it is assumed that the 
relaxation time is so long that the field correlation length 
£z = £,(Atz) = ^{tq/tq) 1 ^ is more or less frozen in until 
the system is again changing slowly, at time At w +Atz- 
The correlation length of the field can only be mea- 
sured indirectly. One of Zurek's proposals, as yet unful- 
filled, is to measure the variance in the superflow in an 
annulus after a quench. Since superflow velocity is pro- 
portional to the gradient of the field phase 9, a random 
walk in phase would suggest that the measurable (AO) 2 
along a perimeter of length L has the form 



(A(9) 2 = O 



L 
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where ^ var measures the effective phase-winding length. 
If, as Zurek does, we assume that £„ ar « £z, then (AO) 2 
is large enough to be observed. 

A more accessible experiment is to measure the density 
of vortices at their formation. If the initial density of 
vortices, the defects oi 4 He, is ndef, with separation £<z e /> 
then 



ndef 



O 
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Zurek makes the assumption that ^def ~ £,z whereby 



ndef 




O -o- 
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Since £o also measures cold vortex thickness, tq 3> t 
corresponds to a measurably large number of widely sep- 
arated vortices. The details do not concern us here, but 
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one experiment § on 4 He, for which t q /t = O(10 10 ), 
shows fair agreement with @, while a second experi- 
ment (as yet unpublished) with comparable tq finds no 
vortices, within its accuracy. Because of this, it is timely 
to reexamine this picture, with its implicit assumption of 
a domain structure characterised by a single length. 

In this brief note we shall argue that the predictions 
(g) and (|]) are unreliable because we cannot identify £ var 
and t;def with £z, even though we may seem to have a 
single-scale theory. To see this we need an explicit model, 
although our conclusions are more general. Firstly, we 
accept that correlation lengths, by whatever definition, 
are frozen in during the time interval — Atz < At < 0, 
and for some time thereafter. In the first instance it is 
then sufficient to perform all our calculations at At = 0. 

At this time the equal-time correlation function is di- 
agonal, 



(0„(r)0 6 (O)) =S ab G(r) 
Suppose that, as Zurek suggested, 



G(r) 



-r/iz 
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for r large compared to £z (but not necessarily asymptot- 
ically large). To go further we assume that, at this time, 
field fluctuations are approximately Gaussian. An imme- 
diate consequence is that the field phases are correlated 
as 



(e i9 ^ e 



m ) « J/(r) 
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for large enough r that f(r) is small, where f(r) = 
G(r)/G(0) is the normalised correlation function, and 
4> a (r) = \4> a (r)\ e* B ( T K Although this seems to support 
a simple domain picture, in which field phases are corre- 
lated on the same scale as fields, the relationship between 
the unobservable £z and the observable £de/ and £ var is 
complicated. 

To begin first with the defect density, the core of every 
vortex is a line zero of the complex field <p. The converse 
is not true since zeroes occur on all scales. However, it 
may be valid Q to count vortices by counting line zeroes 
of an appropriately coarse-grained field, in which struc- 
ture on a scale smaller than £oj the classical vortex size, is 
ignored. That is, we ignore vortices within vortices. This 
is, indeed, the basis of the numerous numerical simula- 
tions of vortex networks built from Gaussian fluc- 
tuations. However, for it to be valid, the result must be 
insensitive to the details of the cutoff. Only then is the 
estimate reliable. For the moment, we put in a cutoff 
I = 0(£o) by hand, as 



G(r) = / fke tk - x G(k)e 
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enabling f(r) to be defined. We stress that if the long- 
distance correlation length £z describes the asymptotic 
behaviour then it depends only on the position of the 



nearest singularity of G{k) in the complex k-plane, inde- 
pendent of I. Even if it only characterises intermediate 
lengths the dependence on I will be weak. 

This is not the case for the line-zero density n zero , 
depending on the short- distance behaviour of G(r), 
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the ratio of fourth to second moments of G{k) e~ k 1 . 
^ Z ero measures the separation of line zeroes. Only if 
dn zero /dl is small at I = £ can we identify n zero with a 
meaningful vortex density, and £ zero with £de/- 

Finally, using the tools of Halperin || and Mazenko et 
al. |TJ] the final correlation length £ var can be shown to 
be given by 
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dr 



i - /M 2 
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where c = 0(1). ^From (|ll| ) we see that the integrand of 
( |l2| ) is O(Cjero) f° r small r, and falls off like exp(—2r/^z) 
for large enough r. If ^ zero = 0(£z) then so is £ va ri but 
otherwise it lies between the two. 

Now that we have established the principle that the 
observable correlation lengths use different attributes of 
the power spectrum G(k) e~ k 1 , we demonstrate it with 
a concrete example, motivated by Zurek's later numer- 
ical |TT| simulations with the time-dependent Landau- 
Ginzburg equation. 

That is, for At w we assume a linear response 



1 d<t> a 

r dt 
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where rj a is Gaussian noise. We go further and assume 
that, for early times, prior to the symmetry breaking at 
least, the self-interaction term can be neglected ((3 = 0). 
This both preserves Gaussian field fluctuations and leads 
to t;z arising in a natural way, as we shall see. [A similar 
approach to relativistic quantum field theory ]10| permits 
a comparison with Kibble's predictions]. All correlation 
lengths depend only on the renormalised /(r), and the 
strength of the noise is immaterial. Writing the resultant 
equation in time and space units tq and £o as 



0a(M) =- 



t 



kf + eo 9(t) 



4>a(k,t) + T r) a (k,t) 
(14) 



gives the solution, at At = 0, </> Q (k) 



r e oTQ 
tq I dtexp 



dt' 



k 2 + [e 6{t) 



%(k). 
(15) 



The resulting un-normalised correlation function has 
power 
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G(k) = / dtexp 



dt' 



-9(f) 
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For a typical quench in 4 He, eo ~ 10~ 2 — 10~ 3 is very 
small, but tq ~ 10 10 is so large that e Tg, e o r Q ^ 1- 
G(fc) can then be approximated as 



G(k) 



T Q ( e0 + fe 2 ) 

' / dte-t'^. 

tq k 2 
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independent of eo. 

On Fourier transforming G(k) of (0) we find 



G(r) oc / dsse 4s e sr cos(sr) 
Jo 
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This does nor give an asymptotic fall-off of the form 
(§). In fact, for large r, G(r) oc exp(-0((r/£z) 4/3 ))- 
Nonetheless, numerically, it is remarkably well repre- 
sented by (|8|), with coefficient unity in the exponent, for r 
being a few multiples of £z > f° r reasons that are not clear 
to us. In that sense Zurek's prediction for a correlation 
length of the form £z = ^(tq/tq) 1 ' 4 is robust. 

However, it is equally easy to determine the density of 
line zeroes at this time. In evaluating 



Szero 



dkk e 



G{k) / J dkk 4 e~ k 1 G{k) 
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we substitute for G{k) from ( |16[ ) and perform the k inte- 
gration first. For small eo, very large tq, and I — O(l), 
in dimensionless units the dominant contribution is from 
t ~ eoTQ. On neglecting terms relatively 0(e~ e ° TQ ) we 
find 



dt 



-t 2 /r Q 



[t + P/2]*/*/ J [t + P/2f/ 



0(1% 
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independent of eo- The details are immaterial. Firstly, 
since £, zer o <C £,z the frozen correlation length of the field 
does not set the scale at which line zeroes appear. If 
To = 10 10 we have O(10 5 ) line zeroes per correlation area. 
Equally importantly, we have a situation in which the 
density of line zeroes depends entirely on the scale at 
which we look. If we look at half the scale (it i.e. half 
a vortex thickness) we see twice as many. It is incorrect 
to identify the length scale £ ze ro with £def since such line 
zeroes cannot be understood as vortices. 

This is not surprising. Although the field correlation 
length £z may have frozen in, the symmetry breaking has 
not been effected. The field magnitude has yet to grow 
to its equilibrium value 



(3/a . 



To see how this happens, we continue to use (13) for 
times t > eoTQ. Then, as the unfreezing occurs, long 
wavelength modes with k 2 < tjrQ — eo grow exponen- 
tially. Provided £qTq ^S> 1 they soon begin to dominate 
the correlation functions. Let 



G n (At) 



dk k n G(k,At) 
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be the moments of G(fc, At), now of the form G(k, At) 

e TQ+At 



dt' exp 

We find that 
G n (At = pAtz) 



reoTQ+At 




- 2 J dt" 


fc 2 + ( 
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dt' 



(23) 



p-(t'-p^TQ) 2 /TQ 

[i' + P/2]™+ 1 /2' 

(24) 



where we measure the time At in units of Atz = ^/tq 

from t = 6qTq and /„ = j Q dkk 2n e~ k . If the linear 
equation ( pi| ) were valid for large p (i.e. large At) then 
the integral is dominated by the saddle-point at t' — 
P^Jtq, to give a separation of line zeroes £ zero (At) of the 
form 



e zero (At = pAt z ) 



G^pAtz) 
G 2 (pAt z ) 



(25) 



independent of the cutoff I. That is, because of the 
transfer of power to long wavelengths, line zeroes become 
widely separated, and £ zero does begin to measure vor- 
tices, and can be identified with £de/- Thus, if the order 
parameter is large enough that it takes a long time (in 
units of Atz) for the field to populate its ground states 
at (|l]) we would recover the Zurek result (^) as an order 
of magnitude result from (p5|), but for entirely different 
reasons. 

Whether we have time enough depends on the self- 
coupling (3, which determines when the linear approx- 
imation fails. For smaller times the integrand gets an 
increasingly large contribution from the ultraviolet cut- 
off dependent lower endpoint, and we recover (|2^). The 
multiple p n of At, at which the exponential modes begin 
to dominate in the moments of G„, can be determined 
by comparing the relative strengths of the contributions 
from the saddlepoint and the endpoint in (24). For the 
former to dominate the latter requires 



Pn 



-1/2 



> 



(rc-1/2) 



-1/2 



(26) 



We see that it takes longer for the long wavelength modes 
to dominate the derivative G2 than the order param- 
eter G\. For tq = O(10 10 ) and I sa 1, pi, satisfying 

3/2 ^ i/4 , 



(21) e Pl /pi > Tq /v27r, must be somewhat larger than 2. 
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Similarly p 2 , satisfying e p ^/p^ 2 > 3^/2tq 4 /^/n is larger 
still, somewhat greater than 4. We stress that these are 
lower bounds. 

On the other hand, at an absolute maximum, the cor- 
relation function must stop its exponential growth at 
At = At sp = qAtz, for some multiple q, when (|^| 2 ), 
proportional to G\, satisfies (|2l|). Let us suppose that 
the effect of the backreaction that stops the growth ini- 
tially freezes in any defects. This then is our prospective 
starting point for identifying and counting vortices. To 
determine q, or At sp , we need to know the strength of 
the noise. In our dimensionless units we have 

( Va (k, t) m (-k',t')) = 6 ab (2k B T){$/a )6(t - t')f{k - k') 

(27) 



(28) 



As a result, (|21j) is satisfied when 

G^qAtz) = Tr 2 a 2 $/f3k B T w 10 2 



for superfluid 4 He at 2K. In order that G\ is dominated 
by the exponentially growing modes, so as to recover 
Zurek's result in the form ( p5| ) , we must have q > P\,P2 
where, from (E8|) 
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3/2 



< 400r, 



1/4 



Q 
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say. In fact, although we need definite figures to make 
comparisons, a coefficient like that in (|29| ) is not entirely 
believable, assuming as it does that the backreaction is 
effectively instantaneous. However, it is probably good 
enough to argue that, for tq = O(10 10 ), q lies between 
Pi and p2- That is, field growth must stop before G 2 
is either well defined, or able to give the result (^). We 
need an impossibly fast quench, with tq = O(10 4 — 10 5 ), 
for q to be large enough. 

What this means for the current He experiments is 
unclear. Although line zeroes will be separated by more 
than a classical vortex thickness they will still show short- 
range structure and are only precursors of vortices. In 
fact, with G\ insensitive to cutoff I, but G2 sensitive, 
£zero is arguably more sensitive to scale than in (19), 
where there was partial cancellation of I dependence be- 
tween numerator and denominator. As for £„ ar the cal- 
culations are much murkier and we have not attempted 
to evaluate it. Much more detailed numerical modelling 
in three dimensions is needed before we can draw any 
quantitative conclusions. In less than three dimensions 



(particularly in one dimension 12 1) the saddle point dom- 
inates more easily over the endpoint, which measures the 
strength of the ultraviolet singularities. 

Even when our approximations are correct, it is not 
clear that the density of vortices at their first produc- 
tion is the relevant extrapolation for experiment. The 
interaction between field modes will tend to redistribute 
the power back to shorter wavelengths in the short-term, 
and the retarded nature of G(r) after At sp will impose 



oscillations that may have consequences, as well as the 
possible incorporation of phonon modes. 

However, we believe that we have presented a con- 
vincing case that, because line-zero density (and phase 
variance) are based on different attributes of the power 
G(k) from the effective long-range correlation length £z, 
we should not make the inference (||) directly. Only if 
there is time enough before the field populates the ground 
states do we recover the Zurek result, although for dif- 
ferent reasons. More details will be given elsewhere [[l3). 
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